Stability and localization of rapid shear in fluid-saturated fault gouge: 1. Linearized stability analysis

Field observations of major earthquake fault zones show that shear deformation is often confined to principal slipping zones that may be of order 1–100 μmwide, located within a broader gouge layer of order 10–100 mm wide. This paper examines the possibility that the extreme strain localization observed may be due to the coupling of shear heating, thermal pressurization, and diffusion. In the absence of a stabilizing mechanism shear deformation in a continuum analysis will collapse to an infinitesimally thin zone. Two possible stabilizing mechanisms, studied in this paper, are rate-strengthening friction and dilatancy. For rate-strengthening friction alone, a linear stability analysis shows that uniform shear of a gouge layer is unstable for perturbations exceeding a critical wavelength. Using this critical wavelength we predict a width for the localized zone as a function of the gouge properties. Taking representative parameters for fault gouge at typical centroidal depths of crustal seismogenic zones, we predict localized zones of order 5–40 μm wide, roughly consistent with field and experimental observations. For dilatancy alone, linearized strain rate perturbations with a sufficiently large wavelength will undergo transient exponential growth before decaying back to uniform shear. The total perturbation strain accumulated during this transient strain rate localization is shown to be largely controlled by a single dimensionless parameter E, which is a measure of the dilatancy of the gouge material due to an increase in strain rate.